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ABSTRACT. A study of the expected frequency of merging galaxies 
is conducted, using the impulsive approximation. Results 
indicate that if we consider mergers involving galaxy pairs 
without halos in a single crossing time or orbital period, the 
expected frequency of mergers is two orders of magnitude below 
the observed value for the present epoch. If we consider mergers 
invovling several orbital periods or crossing times, the expected 
frequency goes up by an order of magnitude. Prelimihary 
calculation indicate that if we consider galaxy mergers between 
pairs with massive halos, the merger is very much hastened. 


INTRODUCTION 

Observations indicate that the frequency of galaxy mergers 
(order of magnitude accuracy only) in the present epoch 
is ~0.3% (Toomre, 1977; Tremaine, 1980), An extrapolation to the 
past yields a frequency ~5%. Loosely bound pairs of galaxies, 
that had separated to great distance in the general cosmic 
expansion, and have lately fallen together again (for the first 
or feweth time) in comet-like plunging ellipses, seem to be the 
most lucrative candidates for mergers. To tackle this problem it 
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is necessary to determine the frequency of merging galaxies on 
the basis of the collision theory and compare the values so 
obtained with the observational ones. 

In this work we have studied sphere-sphere and disk-sphere 
collisions. We have assumed that a galaxy and its nearest 
neighbor are likely to form a bound pair, for frequency 
determinations. We have not embedded the galaxies in massive 
halos, except for one case. A detailed paper on merger of 
galaxies with massive halos will be published soon. The 
determination of general results and the expected frequency of 
merging galaxies is greatly facilitated by using the impulsive 
approximation, as it is a suitable way of studying many 
collisions with varying parameters. As far as the order-of- 
magnitude results are concerned, those obtained by this method 
tally quite well with the results of N-body simulations (Toomre, 
1977; Dekel et al., 1980; Alladin and Narasimhan, 1982; and 
Aguilar and White, 1985). Hence the order of magnitude of the 
frequency that we obtain by this method will be correct, At any 
rate, the crucial questions pertain to the order of magnitude of 
the frequency of merging galaxies. 
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THEORY 


We model the spherical galaxy as a polytrope of index n = 4, 

and the disk galaxy as an exponential model disk with density 

" 4r/R D 2 

distribution, a( r) = a e , where g = M n /( 2irR b) is the 

C v u u 

central density, M D and R D being the mass and radius of the disk 
and b is a constant pertaining to its density distribution. This 
disk is thickened by a method indicated by Rohlfs and 
Kreitschmann (1981). A polytropic bulge of index n = 4 is 
superposed on this disk. The model is discussed in detail in 
Chatter jee, 1989. 

The theory for studying sphere-sphere and disk sphere 
collisions for mergers taking place in a single orbital period is 
basically the same as in Chatter jee, 1987 (hereafter referred to 
as Paper I), except that in the case of disk-sphere collisions 
the appropriate function x = X DS is used, which defines the 
mutual potential energy function for the disk-sphere system (cf. 
Ballabh, 1975). Mj and are the masses and R i = R 2 = R the i r 
common radii. Mergers of disk-sphere systems are dealt with in 
detail in Chatterjee 1990a. 

To study mergers taking place in several orbital periods we 
use the same method as outlined in Chatterjee, 1990b. The method 
is basically a modification of the method used by Alladin, 

1965. In the case of disk-sphere collisions the appropriate 
function x = Xnc used. If r and e denote the polar 
coordinates of the perturber galaxy (Galaxy 1) with respect to 
the test galaxy (Galaxy 2) , in the orbital plane characterizing 
the relative motion, and vfi and E, the reduced mass, the angular 
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momentum and the energy of orbital motion, respectively, then the 


Lagrangian equations for the problem give, 


-1 9 ? 1 / 2 

dt = dr/[2y X {E - W(r ) - r z /( 2 v r z ) } ] 


and de 


-l 2 2 1/2 

Udr)/[2y 1 { E - W(r) - jT/(2yr 2 )}] 


(yr 2 ) 


These equations can be used to determine the orbit and the 
variation of time along the orbit. The relative orbit is studied 
numerically by using these equations. At small intervals of r 
and 0 the tidal effects, and the corresponding change in velocity 
of the stars , is studied by the same method as in Paper I. The 
change in internal energies of the two galaxies AU^ and AU 2 at 
any time t can be obtained from. 


, 1/2 
V(t) = (dr/dt) = [2y 1 {E i - W(t) “ AU^t) - AU 2 (t)}] 


where E. = (1/2) yV. 2 - (GM 1 M 2 /R 2 ) Xgi 


and 

w(t) = w(r ) = vKs) = -(gm x m 2 /r 2 ) x s 
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and are the initial relative velocity and orbital energy of 
the two galaxies, W(t) - W(s) is the mutual potential energy of 
the two galaxies at the instantaneous separation s - 
r/R, x 0 and x c are the values of x = x cc or X n c (depending on 

S ^ S Do L/o 

whether it is a sphere-sphere or disk-sphere collision that is 
under study) (Alladin uses iji = sx) at the initial separation s^ = 
r^/R and s = r/R, respectively. 

We study many collisions between galaxies of equal 
dimensions. In the case of disk-sphere collisions as almost all 
of the mass of the spherical galaxy is concentrated within 1/3 of 
its radius, the effective radius of the spherical galaxy is 1/3 
that of the disk. In the case of mergers taking place in a 
single crossing period, each collision is characterized by a 
value of the distance of closest approach, p, and initial 
velocity V^, as discussed in Paper I, In the case of mergers 
taking place in several orbital periods, the initial separation 
is taken to play the part of p in frequency determination. In 
each case the change in velocity due to dynamical friction is 
calculated for small intervals as the galaxies advance in their 
relative orbit, and the instantaneous relative velocity so 
determined is used to continue the integrations of the equations 
of motion of the two galaxies. Merger takes place when the 
instantaneous relative velocity, V^, of the two galaxies equals 
the velocity of escape between the pair characterized at the 
instantaneous separation, V £ , as discussed in Paper I. In the 
case of mergers taking place in a single crossing time the 
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relative velocity at closest approach, Vp is used for frequency 
determinations, while is taken to play the part of Vp for 
mergers in several orbital periods. 


NUMERICAL RESULTS AND DISCUSSION 


General Resuults 

The computations of many collisions have been carried out, 
each collision being characterized by a value of Vp and p 
(measured in units of the radius of either galaxy) and scaled for 


different mass ratios of the two galaxies. 

For sphere-sphere systems in collision, we find that for 
central impacts, merger takes place in a single crossing time 
provided p < R^y^, the radius containing 1/4 of the mass of the 
victim galaxy, and V <1*2 V R . The merger is affected in a 
single crossing time up to p < R 3 / 4 / the radius containing 3/4 of 
the mass of the victim galaxy; while the merger is effected in 
more than one orbital period up to p ~ R. The merger is effected 
in < 2 orbital periods if R, .. < p < 4R U (4R. * 1/2 of the 

radius for polytropic n = 4 model), R^ being the galactic half- 
mass radius; and in several orbital periods (~ 5 ) for 
4R^ < p < R. In general beyond p ~ 4R^ the dynamical friction 
process becomes quite ineffective and the merger is effected 
extremely slowly. The relative mass of the perturber with 
respect to the victim galaxy does not seem to affect the merger 
times markedly, unless the former has less than about 1/5 the 
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mass of the latter. For p < 4R,, the parabolic perturber orbit 
is circularized in the first orbital period, and the subsequent 
orbit effects the merger; while if p lies between 4R h and R, the 
circularization is achieved in few orbital periods. Beyond 
p ~ R, mergers do not seem to take place as the dynamical 
friction lacks the ability of circularizing the parabolic 
perturber orbit and the same recede to an enormous distance from 
which it does not return in significant orbital timescales. 

In the case of disk-sphere collisions, dynamical friction 
plays a very insignificant role, compared to sphere-sphere 
collisions, due to the small thickness of the victim disk and 
mergers do not seem to be possible for hyperbolic velocities, 
which implies V < V . The merger is effected in a single 
crossing time only if the impact is central or very nearly so, 
and the range of values of p depends strongly on the relative 
mass of the bulge with respect to the disk. For these reasons, 
for the purpose of frequency determinations, it seems to be more 
convenient to determine the range of values of p for which merger 
is possible in single and several orbital periods together (for a 
more detailed treatment see Chatter jee, 1990b). We find that in 
this case mergers are possible up to p ~ R. Only if p > R^ (for 
the victim disk), the merger takes place in’ a much longer time 
than otherwise. Beyond p ~ R the perturber does not seem to 
return in significant orbital timescales. 

We have studied only one collision with galaxies having 
massive halos. So we cannot comment on the frequency of merging 
galaxies with massive halos, but only mention the results of this 
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simulation, as the parameters for this collision have been 
carefully chosen. The model used for the halo has the asymptotic 
dependence M(r) ~ r, where r is the distance measured from the 
center of the galaxy, and is discussed in detail by Allen and 
Martos (1986). This model is excellent for the study of many 
collisions as an analytical treatment of the halo is 
elucidated. The bulge component is modeled as a polytrope of 
index n = 4. A merger of two identical spherical galaxies is 
studied. The halo contains about 10 times the mass of the bulge 
component and 10 its radius, R^. The initial separation of the 
two galaxies is taken to be p = 2R ^ and the initial relative 
velocity is taken to be v^ = 2Vg. Merger takes place in few 
orbital periods (< 3) . The exact number of orbital periods 
depends upon the value of the separation of the centers of the 
two galaxies for which merger is defined to be complete. For 
frequency determinations, whether or not the merger takes place 
in a significant orbital timescale is the crucial fact. 

Brief Comparison with Previous Work 

The results of van Albada and van Gorkom (1977), White 

(1978), Roos and Norman (1979), are summed up in Figure 1 of 

Aarseth and Fall (1980). Translated in terms of our parameters, 

this figure is indicative of a high frequency of mergers in the 

region corresonding to p < 0.8 x 2R^ = 0.216R, for polytrope n = 

4 model, and V_ < V < 1.2V_. Villumsen (1982) finds that for the 
E ~ p ~ E 

distance of closest approach p = R^ < 10 units, corresponding to 
0.22 R (as the tidal radius for his model is 45 units), merger is 
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quite frequent; while for p = = 10 units, merger occurs only 

for low relative velocities. Alladin and Parthasaraty (1978) 
indicate that stellar systems, modeled as polytropes of index 
n = 4, in mutual circular orbit, merge in about 15% of an orbital 
period for p < 0.2R; while for p > 0.2R, merger takes 
considerably more time. All these results compare well with our 
results that merger in a single crossing time ceases, for sphere- 
sphere pairs, for p < ~ 0.2R (for polytrope n = 4). The 

results are not very much model dependent. 

In the region R 3/4 < p < R^, merger is achieved in a 
time ~ 2 orbital periods, while it is achieved in several orbital 
periods in the region 4R^ < p < R. These results compare quite 
well with the results of Borne (1984). (For more details see 
Chatterjee, 1990b.) It is interesting to note that Borne 
includes a relaxation time in his scheme, which we do not, but 
the resonating stars which hasten the orbital decay of the 
galaxies in his simulations play a comparatively insignificant 
part in our scheme. 

Frequency Determinations 

In dense regions a galaxy and its nearest neighbor can be 
visualized to form a loosely bound pair (Toomre, 1977; Tremaine, 
1980). As mentioned before, the general opinion is that 
collisions between initially loosely bound pairs lead to the 
majority of mergers. In dense regions the average distance 
between galaxies is ~ 10R a , where R a is the average radius of a 
galaxy (Mitton, 1977; Ogorodnikov, 1965). So though 
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its mean 


theoretically the parameter p can vary from 0 to », 
value in regions where mergers are frequency, can be taken 
as p <10 (in units of the radius of a galaxy). We take the 

OO 

total range of variation of V p to be from 0 to 4V E , as collisions 

with velocity V > 4V„ seldom occur 
J p ~ E 

(i) Mergers of Spherical Systems in a Single Crossing Time. 

Central Impacts: p < <1*2 V £ . The probability in terms 

of the energy of the collision is given by 

P„ = (1.2 V„ ) 3 /( 4V,., ) 3 = 0.027; and the probability in terms of 

h hi ej 

2 -5 

the impact parameter is given by P^ = ( R ]y 4 P w ) = 4.56 x 10 ; 

which gives the probability of mergers in this case as P c = P E x 
P = 1.2 x 10” 6 . Which gives the expected frequency of such 

Jr 

mergers as 1.2 x 10 “^ %. 

Off Center Impacts: R. < p leading to mergers only 

1/4 ~ ~ o/4 

if Vp < V £ . The corresponding probabilities in this case are 
given by P E = [V E /(4V E )] 3 = 0.015625, P p = (R 2 3 / 4 - (R 2 i/ 4 )/p 2 = 
3.645 x 10“^, giving P oc = P E x P p = 5.7 x 10”^. This gives an 
expected frequency of 5.7 x 10“^ %. Since these are mutually 
exclusive events, the probability of mergers of spherical systems 
in an orbital period is, P = P c + P QC = 6.9 x 10”^, which leads 
to an expected frequency of such mergers as 6.9 x 10“^ %. 

(ii) Mergers of Spherical Systems in Several Orbital Periods. 
R 3/4 < p < R, V <V E . In this case the respective 

probabilities are given by P E * [V E /( 4V E __] 3 = 0.015625, P p = 
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( r2 3/4 “ r2 i/ 4)/P 2 = 9 « 6 x 10-3, which gives P = P E x P p = 1.5 x 
10”^; which gives the expected frequency of such mergers as 1.5 x 
10" 2 %. 

Hence the expected frequency of mergers of spherical systems is 
given by 1.57 x 10“ 2 %. 

(iii) Mergers of Disk-Sphere Systems. 

p < R, Vp < Vg. In this case the respective probabilities are 
given by P E [V E /(4V E )]3 = 0.015625, P p = (R/p^) 2 = 0.01 , which 
gives the merger probability as P = P E x P p = 1.56 x 10“^, which 
gives the expected frequency of such mergers as 1.56 x 10“ 2 
% - 10 -2 %. 

CONCLUSIONS 

In general we conclude that the expected frequency of 
merging galaxies (not considering the effect of massive halos) 
is ~ 0.01%, which is an order of magnitude below the 
observational value ~0.1 %. The expected frequency of merging 
galaxies increases by an order of magnitude when we take into 
account mergers occurring in several orbital periods, as compared 
to those occurring in a single crossing time. More than 90% of 
the mergers take place in several orbital periods, and there is 
an indication that the majority of them take place in a 
time ~2 to 3 orbital periods. For mergers in a single crossing 
time, more than 80% of them are due to off-center impacts. 
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Spherical and disk-sphere systems in merger seem to be equally 
frequent. 

If our single simulation with galaxies embedded in massive 
halos is at all typical (which in our view it is), then the 
frequency of galaxy mergers with massive halos will be at least 
an order of magnitude higher than in the absence of halos 
(approaching the observational value for the present epoch). We 
are studying more collisions of type to throw light on this 
aspect . 
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